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As originally formulated, the Generalized Alignment Index (GALI) method of chaos detection 
has so far been applied to distinguish quasiperiodic from chaotic motion in conservative nonlinear 
dynamical systems. In this paper we extend its realm of applicability by using it to investigate 
the local dynamics of periodic orbits. We show theoretically and verify numerically that for stable 
periodic orbits the GALIs tend to zero following particular power laws for Hamiltonian flows, 
while they fluctuate around non-zero values for symplectic maps. By comparison, the GALIs of 
unstable periodic orbits tend exponentially to zero, both for flows and maps. We also apply the 
GALIs for investigating the dynamics in the neighborhood of periodic orbits, and show that for 
chaotic solutions influenced by the homoclinic tangle of unstable periodic orbits, the GALIs can 
exhibit a remarkable oscillatory behavior during which their amplitudes change by many orders 
of magnitude. Finally, we use the GALI method to elucidate further the connection between 
the dynamics of Hamiltonian flows and symplectic maps. In particular, we show that, using for 
the computation of GALIs the components of deviation vectors orthogonal to the direction of 
motion, the indices of stable periodic orbits behave for flows as they do for maps. 
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1. Introduction 

The method of the Generahzed Ahgnment Indices (GALIs) was originahy introduced in |Skokos et a/.L 
20071 ] as an efficient chaos detection method. To date, the GALI method has been successfuhy apphed 



to a wide variety of conservative dynamical systems for the discrimination be tween regular and chaotic 
motion, as well as for the detection of regular motion on low dimensional tor i lAntonopou los fc Bountis 



20061; IChristodouhdi fc Bountisl. [20oi; ISkokos et^, l2QQ8|; [Manos et al 



Skokos fc Gerlachl . l2010l : iManos fc RAiffol . l2010l : iManos fc Athanassoulal . 



2008a 



2011 



2009; Bountis et a/J . 120091 
Gerlach et a/] . l201l[ ]. 



In the present paper we extend and complete the study of the GALI method, focusing on its behavior for 
the special case of periodic orbits and their neighborhood in conservative dynamical systems. The detection 
of periodic orbits and the determination of their stability are fundamental approaches for the study of 
nonlinear dynamical systems, since they provide valuable information for the structure of their phase 
space. In particular, stable periodic orbits are associated with regular motion, since they are surrounded 
by tori of quasiperiodic motion, while in the vicinity of unstable periodic orbits chaotic motion occurs. 

The GALI method is related to the evolution of several deviation vectors from the studied orbit, and 
therefore is influenced by the characteristics of the system's tangent space. The main goal of the paper is 
to determine the usefulness of the method for probing the local dynamics of periodic orbits with different 
stability types. We manage to achieve this goal by deriving theoretical predictions for the behavior of 
GALIs for stable and unstable periodic orbits. We also verify numerically the validity of these predictions, 
by studying the evolution of GALIs for periodic orbits of several Hamiltonian flows and symplectic maps, 
clarifying also the connections of such dynamical systems. In addition, we show how the properties of 
the index can be used to locate stable periodic orbits, and to understand the dynamics in the vicinity of 
unstable ones. 

The paper is organized as follows: in the first two introductory sections we recall the definition of the 
GALI, describing also its behavior for regular and chaotic orbits (Sect. [2]), and report the several stability 
types of periodic orbits in conservative systems (Sect.[3j). In Sect.[4]we first study theoretically the behavior 
of the index for stable and unstable orbits, and then present applications of the GALI to particular orbits 
of Hamiltonian flows and symplectic maps. Sect. [5] is devoted to the dynamics in the neighborhood of 
periodic orbits, while Sect. [6] is dedicated to the relation between the GALIs of stable periodic orbits for 
flows and maps. Finally, in Sect. El we summarize our results. 



2. The Generalized Alignment Index (GALI) 

Let us briefly recall the definition of the GALIs and their behavior for regular and chaotic motion in 
conservative dynamical systems. Consider an autonomous Hamiltonian system of N degrees of freedom 
(A^dof), described by the Hamiltonian H{qi, q2, . . . , qNiPi,P2, • • • iPn), where Qi and p^, i = 1, 2, . . . , are 
the generalized coordinates and conjugate momenta respectively. An orbit in the 2A^-dimensional phase 
space S of this system is defined by a vector x{t) = {Qi{t),q2{t), • • - ^qN{t)^Pi{t),P2{t), . . . ,pN{t))^ with 
Xi — qi^ Xi^N = Pi, ^ = 1, 2, . . . , A^. The time evolution of this orbit is governed by Hamilton's equations of 
motion 

dx ^^H dH\ 



while the time evolution of an initial deviation vector w(0) = ((ixi(O), . . . , dx2N{0)) from the x{t) solution 
of Eqs. ([1]), obeys the variational equations 

^ = M(f(t)).^, (2) 

where M = dV /dx is the Jacobian matrix of V. 

Let us also consider a discrete time t = n G N conservative dynamical system defined by a 2A^- 
dimensional (2A^D) symplectic map F. The evolution of an orbit in the 2A^-dimensional space S of the 
map is governed by the difference equation 

x{n + 1) = Xn+l F{xn)- (3) 
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In this case, the evolution of a deviation vector w{n) 
by the corresponding tangent map 

w{n + 1) = Wn+l 



with respect to a reference orbit x^, is given 



(4) 



For A^dof Hamiltonian flows and 2A^D maps the Generahzed Ahgnment Index of order k (GALI/^), 
2 < k < 2N, is determined through the evolution of k initially linearly independent deviation vectors 
Wk(0). To avoid overflow problems, the resulting de viation vector s i^J/^ft ) are continually normalized, but 
their directions are kept intact. Then, according to [S kokos et al . [2007l | GALI/^ is deflned as the volume 
of the fc-parallelogram having as edges the k unit deviation vectors Wi{t) = Wi(t)/\\wi(t)\\, i = 1, 2, . . . , fc, 
determined through the wedge product of these vectors as 

(5) 



GALIfc(t) = \\wi{t) ^Mt) ^Mt)l 



with II • II denoting the usual norm. From this deflnition it is evident that if at least two of the deviation 
vectors become linearly dependent, the wedge product in Eq. ([5]) becomes zero and the GALI/^ vanishes. 

In the 2A^-dimensional phase space S of an A^dof Hamiltonian flow or a 2A^D map, regular orbits lie 
on 5-dimensional tori, with 2 < s < N for Hamiltonian flows, and 1 < 5 < for maps. For such orbits, 
all deviation vectors tend to fall on the 5-dimensional tangent space of the torus on which the motion 
lies. Thus, if we start with k < s general deviation vectors, these will remain linearly independent on the 
5-dimensional tangent space of the torus, since there is no particular reason for them to become linearly 
dependent. As a consequence GALI/^ remains practically constant and different from zero for k < s. On the 
other hand, GALIj^ tends to zero for /c > 5, since some deviation vectors will eventually have to become 
linearly depende nt. In particular, the generic behavior of GALlk for regular orbits lying on 5-dimensional 



tori is given by [ Christodoulidi fc Bountisl . l2006l : ISkokos et all 120081 ] 



GALIfc(t) oc 



constant if 2 < A: < 5 
^ ifs<k<2N -s 
1 if2N -s<k<2N 



(6) 



Note that these estimations are valid only when the conditions stated above are exactly satisfled. For 
example, in the case of 2D maps, where the only possible torus is an 5 = 1-dimensional invariant curve, 
the tangent space is 1-dimensional. Thus, the behavior of GALI2 (which is the only possible index in this 
case) is given by the third branch of Eq. ([6]), i.e. GALI2 oc since the first two cases of Eq. ([6j) are not 
applicable. From Eq. (|6]) we deduce that, for the usual case of regular orbits lying on an A^-dimensional 
torus, the behavior of GALI^^^ is given by 



^ , ^ ^ , , f constant if 2 < k < N 
GALIfc(t) oc <^ 1 - - 



if N <k<2N 



(7) 



On the other hand, for a chaotic orbit all deviation vectors tend to become linearly dependent, aligning 
themselves in the direction defined by the maximum Lyapunov char acteristic exponent (mLCE) and hence, 
in that case, GALI^^ tends to zero exponentially following the law [Sko kos et all |2007[ | 

GALIfc(t) oc e-[(^i-^2)+(^i-^3)+-+(ai-^.)]t^ (8) 

where ai, . . . , cr/^ are the first k largest Lyapunov characteristic exponents (LCEs) of the orbit. 

The GALI is a generalization of a similar indicator called the Smaller Alignment Index (SALI) 
Skokosl . l2001bl : ISkokos et all 12003 ^ 2004 ] , which has been used succes s fully f or the detection o f chaos 



in several dynamical systen i s ISzell et all 12004 



PanagODOulos et al . 2004 : Bountis fc Skokosl, [2OO6 



Capuzzo-Dolcetta et all l2007l : iManos et all l2008bl : iMacek et all l2007l : IStranskv et all l2009l : iMacek et al 



20101]. The generalization consists in the fact that the GALIs use information of more than two deviation 



vectors from the reference orbit, leading to a faster and clearer distinction between regular and chaotic 
moti on compared with S ALI. In practice, SALI is equivalent to GALI2 since GALI2 oc SALI (see Appendix 



B of jSkokos et all l2007l | for more details). 



For the numerical computation of GALIs we consider the k x 2A^ matrix W(t) having as rows the 
coordinates Wij{t) of the unit deviation vectors Wi{t), z = 1, 2, . . . , A:, j = 1,2, ...,2A^, with respect to 
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the usual orthonormal basis ei = (1, 0, 0, ... , 0), 62 = (0, 1, 0, ... , 0),..., e2N — (0, 0, 0, ... , 1) of the 2N- 
dimensional phase space S. Thus, GAlAj^{t) can be evaluate d as the square roo t of the sum of the squares 
of the determinants of all possible k x k submatrices of W [Skokos et al. . 20071 ] 



GALK 



E 

l<ii<i2<---<ifc<2A^ 



2\ 



1/2 



(9) 



Here the sum is performed over all possible combinations of k indices out of 2N, \ • | denotes the determinant, 
and the explicit dependence of all quantities on the time t is omitted for simplicity. 

Equation ([9]) is ideal for the theoret ical determination of the a s ymptotic behavior of GALIs for chaotic 
and regular orbits. It has been used in jChristodoulidi fc Bountisl . I2OO6I : ISkokos et all I200I I2OO8I ] for the 
derivation of Eqs. ([6j) and ([8]), and will be applied later on in Sect. 14.11 for the determination of GALIs' 
behavior for periodic orbits. However, from a practical point of view the application of Eq. ([9]) for the 
numerical evaluat ion of GALI/^ is not very efficien t as it might require the computation of a large number of 
determinants. In [ Antonopoulos fc Bountis . I2OO6I : ISkokos et aZI . I2OO8I ] a more efficient numerical technique 
for the computation of GALI^, which is based on the Singular Value Decomposition of matrix W(f) was 
presented. In particular, it has been shown that GALI^ is equal to the product of the singular values Zi > 0, 
i = 1,2,... ,fc of W^{t) 



GALh{t) = l[zi{t), 



(10) 



where (^) denotes the transpose matrix. 



3. Stability of Periodic Orbits 

Now, consider a T-periodic orbit (i.e. an orbit satisfying x{t -\- T) = x{t)) of an A^dof Hamiltonian flow 
or of a 2A^D symplectic map. Its linear stability is determined by the eigenvalues of the so-called mon- 
odromy matrix Y(T), which is obtained from the solution of th e variation al equations for one period T 
(see for examp le |Broucke ^ 19691 ] Lichtenberg fc Liebermanl . ll992L Sect. 3.3] [S kokoa . 12001 al : iHadjidemetrioul . 
l2006[ | jCvitanovic et a/ll2009l . Chapt. 4, 5]). The monodromy matrix is symplecticG, and its columns corre- 
spond to linearly independent solutions of the equations that govern the evolution of deviation vectors. In 
particular, the evolution of an initial deviation w{0) from a T-periodic orbit is given by 



wiiT)^[Y{T)Y-w{0), 



1,2, 



(11) 



Due to the symplectic nature of the monodromy matrix and the fact that its elements are real, the 
eigenvalues of Y(T) have the following property: if A is an eigenvalue then 1/A and the complex conjugate 
A* are also eigenvalues. This property shows that the eigenvalues A = 1 and A = — 1 come in pairs and that 
complex eigenvalues with modulus not equal to 1 always appear in quartets. When all eigenvalues are on 
the unit circle the corresponding periodic orbit is said to be stable. If there exist eigenvalues off the unit 
circle the periodic orbit is unstable. 

A few remarks on the connection of Hamiltonian systems with symplectic maps are necessary at this 
point. Since autonomous Hamiltonian systems are conservative, the constancy of the Hamiltonian function 
introduces a constraint which flxes an eigenvalue of the monodromy matrix to be equal to 1 and so, by 
the symplectic property, there must be a second eigenvalue equal to 1. Thus, for an A^dof Hamiltonian 
system there are only 2(A^ — 1) a priori unknown eigenvalues, and so we can reduce our study to a 2(A^ — 1)- 
dimensional subspace of phase space S. This subspace is obtained by the well-known method of the Poincare 



^Y(T) satisfies the condition Y(T)^ • J2N • Y(T) 
and Otv is the N x N zero matrix. 



J2N1 with J 



2N 



, where I^v is the N x N identity matrix 
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surface of section (PSS) (e.g. Lichtenberg fc Liebermanl . Il992l . p. 17-20] [Cvitanovic et al . I2OO9I . Sect. 3.1, 
3.2]). The corresponding monodromy matrix of the periodic orbit is also symplectic. Thus, in this sense, 
an A^dof Hamiltonian system is dynamically equivalent to 2(A^ — 1)D symplectic map. 

The different stability types of a periodic orbit in H amiltonian systems of 2dof and 3dof (or equivalent! 
in 2p and 4D maps) have b een studied in detail in [Br ouckdJ I969l : 1lladj idemetrioul . [l975l : IPullin fc Meisi 



19981 : Hadj idemetrioul . 2006], while the stabi l ity of periodic orbits in hig:her dimens i onal co nservative sys- 
tems was considered in [Howard fc MacKayl . Il987l : iHoward fc Dullinl . Il998l : ISkokod . 12001^ . Following the 
terminology introduced in [S koko'sn2001al ]. the general stability type of a periodic orbit of an A^dof Hamil- 
tonian system, or equivalently a 2(A^ — 1)D map, is denoted by 



SpUjj 



Ah with p + m + 2/ = 7V-l, 



(12) 



which means that p couples of eigenvalues are on the unit circle, m couples are on the real axis and / 
quartets are on the complex plane but off the unit circle and the real axis. We conclude that a periodic 
orbit is stable only when its stability type is Sn-i- In all other cases the orbit is unstable since there exist 
eigenvalues of the monodromy matrix off the unit circle. For example, in the case of a 3dof Hamiltonian 
system or a 4D map a periodic orbit can be linearly stable (^2) or have three different types of instabil- 
ity: S^U^, U2, Ai (often called simp le unstable, double unstable and complex unstable respectively, see 



e.g. ( Contopoulos fc Magnenatl . 1985l |). 



4. The Behavior of the GALI for Periodic Orbits 
4.1. Theoretical treatment 



Let A^, i — 1,2,... ,2A^ be the (possibly complex) eigenvalues of the monodromy matrix Y(r) of a T- 
periodic orbit, ordered as |Ai | > IA2 I > ••• > |A2a^|. Then, the corresponding LCEs (jj, i = 1 , 2, . . . , 2N are 



(13) 



given by jBenettin Galganij . 119791 : iBenettin etall . 119791 : ISkokos et oR l2007l : ISkokosl . l2010l ] 



a, 



In the case of unstable periodic orbits, where at least |Ai| > 1, we get ai > 0, which implies that 
nearby orbits diverge exponentially from the periodic trajectory. Unstable periodic orbits of non-integrable 
Hamiltonian systems and symplectic maps are located inside chaotic domains. All non-periodic chaotic 
orbits in these domains have the same spectrum of LCEs, which in general differs from the spectrum of 
LCEs of the unstable periodic orbits of these domains. 

For de termining the behavior of GALIj^^ for unstable periodic orbits, one can apply the analysis pre- 
sented in [Skoko s et all 12007 ] for chaotic orbits which also have ai > 0. This approach leads to the 
conclusion that GALIj^^ of unstable periodic orbits tends to zero exponentially following the law ([8]) 



GAUkit) oc e-[(^i-^2)+(c7i-a3)+-+(ai-afe) 



(14) 



However, the case of stable periodic orbits needs a more careful investigation. For this purpose let 
us consider an A^dof Hamiltonian system expressed in action-angle variables J^, i = 1,2, ... ,A^. The 
equations of motion of a periodic orbit of this system are 

dH ^ • dH 



0, 



^i(^i,^2,...,^iv), 1 < i < A^. 



d0i ' ' dJi 
The frequencies cji satisfy a relation of the form 

ki k2 " ' kN 

where A:^, z = 1, 2, . . . , A^, are integer numbers and ri( Ji, J25 • • • , -^A^) 
the orbit. Equations (fTSll can be easily integrated to give 

Ji{t) = Jio, 9i{t) = 9io + n{Jio,J2o,---,JNo)kit, l<i<N 

where J^o, ^io? ^ = I5 2, . . . , A^ are the initial conditions. 



(15) 



(16) 



27V /T with T being the period of 

(17) 
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Let us now denote by 77^, i — 1,2, ... ,A^, small deviations from Ji and Oi respectively. Inserting 
Eqs. (ITSjl and (1T6]) into the variational equations of the Hamiltonian system we get 



N 



= 0, f],^hY,^3^3^ l<i<N, 



(18) 



where — dVt/dJj are computed for the initial constant values Jjo, j = 1, 2, . . . , A^. Using as basis of the 
2A/^-dimensional tangent space of the Hamiltonian flow the 2N unit vectors {^i, 1)2, • • • , V2n}^ such that the 
first N of them correspond to the N action variables and the remaining ones to the N conjugate angle 
variables, any initial deviation vector Wi{Q) — (CU^)' ^2(0)' • • • ' CVl^)? '^i(O)' ^2(0)' • • • ' ^7v(0))' evolves in 
time as 



N 



N 



N 



(19) 



From the above it readily follows that for sufficiently long times ||i?i(t)|| oc t. 

Let us now consider fc, initially linearly independent, randomly chosen, unit deviation vectors 
{wi, . . . , Wk}i with 2 < k < 2N, and let W be the matrix having as rows the coordinates of these vectors 
with respect to the {vi,V2, • • • , V2n} basis. Defining by and ry^, z = 1, 2, . . . , A^, the A: x 1 column matrices 



Vi 



(20) 



W(t) assumes the form 




N 



kit 



N 



k]s[t 



(21) 



W(t)oci.W'^(t) = l 

where we have considered Y\!i= 

In order to determine the leading order behavior of GALI/^ as t grows, we look for the fastest increasing 
determinants of all A: x /c minors of matrix W^. For 2 < k < 2N — 1, these determinants include only one 
column of containing the term ^l^i 



Wi 



(X t''. Then, Eq. © can be used for the computation of GALIfc. 



than one columns proportional to 
t-(fc-i) for 2 < A; < 2Ar-l. For k = 



and grow proportional to since determinants with more 

are identically zero. Thus, we conclude that GALlk(t) oc 

2N, is a square 2N x 2N matrix which has a constant determinant, 
since time appears only through multiphcations with the N first columns of W^, and so GALIj[^(t) (x 
Summarizing, the time evolution of GALI/^ for stable periodic orbits of A^dof Hamiltonian systems is given 

by 

if 2 < < 2iV - 1 



GALIfc oc 



iik = 2N 



(22) 



It is worth mentioning that Eq. (l22l) can be retrieved from Eq. (|6]) by assuming motion on an 5 = 1- 
dimensional torus, i.e. on an 1-dimensional curve, which is the stable periodic orbit. Note that for 5 = 1, 
only the last two branches of Eq. (|6j) are meaningful. 

Stable periodic orbits of symplectic maps correspond to stable fixed points of the map, which are 
located inside islands of stability. Any deviation vector from the stable periodic orbit performs a rotation 
around the fixed point. This, for example, can be easily seen in the case of 2D maps where the islands in 
the vicinity of a stable fixed poin t can be re presented through l ineari zation, by ellipses (see for instance 
Lichtenberg fc Liebermanl . Il992l . Sect. 3.3.b] |Lega fc Froeschlel . \200l\\ ). Thus, any set of 2 < /c < 2A^ 



initially linearly independent, unit deviation vectors will rotate around the fixed point, keeping on the 
average the angles between them constant. This means that the volume of the fc-parallelogram having as 
edges these vectors will remain practically constant, exhibiting some fluctuations, since the rotation angles 
are constant only on average. So, in the case of stable periodic orbits of 2A^D maps we have 



GALIfc (X constant, 2 <k < 2N. 



(23) 
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4.2. Numerical results - Hamiltonian flows 

To verify the validity of the theoretical predictions of Eqs. (fHl) and ([22D we now compute the GALIs for 
some representative Hamiltonian systems of different number of degrees of freedom. 



4.2.1. 2dof Henon-Heiles system 

First we consider the well-known 2dof Henon-Heiles model 



Henon fc Heilesl . 119641 ] 



H2 = + pD + ^(^' + y") + ^"y - ^y"' (24) 

In our study we keep the value of the Hamiltonian fixed at H2 — 0.125. Fig. [U^a) shows the PSS of the 
system defined by x = 0, > 0. We consider two stable periodic orbits (whose stability type is Si 
according to Eq. (fT2ll ): An orbit of period 5 (i.e. an orbit intersecting the PSS at the 5 points denoted 
by blue crosses in Fig. [T](a)) with initial condition (x,y,Px,Py) ^ (0.0,0.35207,0.36427,0.14979), and an 
orbit of period 7 (red squares in Fig. [11(a)) with initial condition (x, 

y 5 Px 1 Py ) 

(0.0,0.45882,0.32229,0.0). 

The time evolution of GALI/^, k = 2,3,4: for these two orbits, for a random choice of initial orthonormal 
deviation vectors, is shown in Figs. [11(b) and (c) respectively. For both orbits the indices show a power law 
decay to zero, in accordance with the theoretical prediction of Eq. (|22l) for N — 2. 





100 200 300 400 500 600 1 2 3 4 5 6 

t Log(0 



Fig. 1. (a) The PSS of the 2dof Henon-Heiles system ()24p with H2 — 0.125. The intersection points of stable periodic orbits 
of period 5 (blue crosses) and 7 (red squares), as well as an unstable orbit of period 5 (green circles) are also plotted. The 
line Py — 0^ denoting a set of initial conditions discussed in Sect.[5l is also plotted. The time evolution of GALI2 (red curves), 
GALI3 (green curves) and GALI4 (blue curves) for these three orbits is shown in panels (b), (c) and (d) respectively. Both 
axes of (b) and (c), and the vertical axis of (d) are logarithmic, (e) The time evolution of the quantity Li(t), which has as 
limit for t ^ 00 the mLCE ai of the unstable periodic orbit (horizontal dotted line). Plotted lines correspond to functions 
proportional to t~^ ^ and t~^ in (b) and (c), and the exponential laws ([14]) for ai = 0.084, (72 = in (d). 



In order to check the validity of Eq. (fT4l) , we consider an unstable periodic orbit (of 
Ui type) of period 5 (green circles in Fig. [11(a)) with initial condition {x,y,px^Py) ^ 
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(0.0,0.2083772012,0.4453146996,0.1196065752). The theoretically expected value of this orbit's mLCE 
(Ji is estimated from Eq. (|T3l) to be ai ^ 0.084, while ^2 = because the Hamiltonian function is an 
integral of motion. 

In Fig. [T](d) the time evolution of the corresponding GALI^^^, = 2,3,4 is plotted. From these results 
we conclude that the computed values of GALIs are well approximated by Eq. (fT4|l for ai = 0.084 and 
(72 = 0, at least up to t ^ 350. After that time we observe a change in the exponential decay of GALI2. 
This happens because the numerically computed orbit deviates from the unstable periodic orbit, due to 
computational inaccuracies, and enters the surrounding chaotic domain, which is characterized by different 
LCEs. This behavior is also evident from the evolution of the finite time mLCE Li{t) (Fig. [U^e)) having as 
limit for t ^ oc the mL CE ai o f the c omputed orbit (for more details on the computation of the mLCE 
the reader is referred to Skokosl . 2010l . Sect. 5]). For an initial time interval, Li{t) well approximates the 



mLCE of the unstable periodic orbit, but later on, due to the divergence of the computed orbit from the 
periodic trajectory, Li{t) tends to a different value, which is the mLCE of the chaotic domain around the 
unstable periodic orbit. 

4.2.2. A 3 do f Hamiltonian system 

Let us now investigate the behavior of the GALIs for a 3dof Hamiltonian system, where different types of 
unstable periodic orbits can appear. In particular, we consider a system of three harmonic oscillators with 
nonlinear coupling, described by the Hamiltonian 

The harmonic frequencies of the oscillators are determined by parameters A, C, and the strengths 
of the nonlinear couplings by e and 77. This system was introduced as a crude description of the inner 
parts of distorted 3-dimensional elliptic g:alaxies. Detailed studies of its basic families of periodic orb its 
were performed in |Contopoulos fc Barbanisl . 1 19851 : IContopoulos fc Magnenatl . 1 19851 : IContopoulosl . Il9868 



Hz = i^ipl+pl+vl) + i^{Ax^ + By'' + Cy'') - exz' - r^yz\ (25) 



Following these works, we fix A = 0.9, B = 0.4, C = 0.225 and = 0.00765 and vary e and 77 in order to 
study periodic orbits of different stability types. 

In Fig. [2]^a) we plot the time evolution of GALIs for a stable {S2) periodic orbit with initial condition 
{x,y,z,px,Py,Pz) ^ (-0.06686,0.01230,0,0,0,0.10590) for e = 0.2 and r] = 0.1. The 3dof system has a 
6-dimensional phase space and so, 5 different GALIj^, with 2 < A: < 6, are defined. All GALIs decay to zero 
following the power law predictions given by Eq. (122]) for = 3. 

Let us now study representative cases of all the different types of unstable periodic orbits that can 
appear in a general 3dof system. In particular, we consider an SiUi periodic orbit with initial condition 
{x,y,z,p^,Py,p,) ^ (-0.0238841214,0744533850,0,0,0,0.1121127613) for e = 0.848, 7/ = 0.1 (Fig. [2(b)), 
an U2 periodic orbit with initial con- 

dition {x,y,z,p^,py,p,) ^ (-0.0392937629,0.0648373644,0, -0.0564496390,0.0021636015,0.0950663122) 
for £ = 0.35, 77 = 0.51 (Fig. El^c)), and a Ai periodic orbit with initial condition {x,y,z,Px,Py,Pz) ^ 
(-0.0456720106,0.0658047594,0,0,0,0.1081228661) for e = 0.6 and r] = 0.3 (Figs. [2(d) and (e)). 

Using Eq. I^Eij we estimated the LCEs to be ai ^ 0.046, cr2 = and ai ^ 0.014, (72 ^ 0.0019 for the 
SiUi and the U2 unstable periodic orbits respectively. Using these values as good approximations of the 
actual LCEs, we see in Figs. [2(b) and (c) that the evolution of GALIs is well reproduced by Eq. ([T4|) . 

An eigenvalue of the monodromy matrix of the Ai unstable periodic orbit is numerically found to be 
Al ^ 1.410 + 0.164z, while the remaining three of them (apart from the two unit ones) are 1/Ai, and 
1/A^. Then, from Eq. (|T3]1 we estimated the three largest LCEs of the periodic orbit to be ai = cr2 ^ 0.023, 
(73 = 0. The evolution of the GALIs for this orbit is shown in Fig. [2(d). Although the periodic orbit is 
unstable, GALI2 does not decay to zero but remains constant until t ^ 10^. This happens because, according 
to Eq. ([M]) GALI2 oc e-(^i-^2)t^ but in this case o"! = cr2. However, due to unavoidable inaccuracies in the 
numerical integration, the computed orbit eventually diverges from the unstable periodic one and enters 
a chaotic domain characterized by different LCEs with ai 7^ (72. This divergence is also evident from the 
evolution of q uantities Li(t), L2{t) in Fig. [2(e), whose limits at t ^ 00 are ai and (72 respectively (see 
Skokosl . [20T0I for more details on the computation of ai and (72). In particular, we get Li(t) ^ L2{t) for 
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-{2 <x,-c72)/lnlO 
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Log(0 





Fig. 2. The time evolution of GALI/., 2 < /c < 6 for an (a) S2 stable, (b) SiUi unstable, (c) U2 unstable and (d) Ai 
unstable periodic orbit of the 3dof Hamiltonian system ([25)) . Both axes of (a) and (d), and the vertical axis of (b) and (c) are 
logarithmic. Plotted lines correspond to appropriate power laws {|22]) in (a), and exponential laws (|14|) in (b) and (c). (e) The 
time evolution of quantities ivi(t), L2{t) having respectively, as limit the two largest LCEs ai, cf2 of the Ai unstable periodic 
orbit. The theoretically estimated value ai = (72 = 0.023 is denoted by a horizontal line. 



t < 10^, while later on the two quantities attain different values. Consequently, for t > 10^ GALI2 starts to 
decay exponentially to zero. On the other hand, all other GALIs in Fig. [2(d) show an exponential decay, 
even when GALI2 remains constant, since the corresponding exponents in Eq. (IT4|) do not vanish. 



4.2.3. A multi- dimensional Hamiltonian system 

Finally, we turn to a multi-dimensional Hamiltonian system representing a 1-di mensional cha i n of 5 identical 
particles with nearest neighbor interactions given by the FPU-/? Hamiltonian Fermi et al . 1955[ | 



1 ^ ^ / 1 1 



j=l j=o ^ ^ 

where xj is the displacement of the jth particle from its equilibrium position and pj is the corresponding 
conjugate momentum. In our study, we set /3 = 1.04 and impose fixed boundary conditions to the system, 
so that we always have xq = xq = 0. 

Let us consider t wo pa rticular members of a family of periodic orbits studied in [Ooyama et al I, Il969l : 
lAntonopoulos et al. I 120061 1 which have initial conditions of the form xi(0) = —xs{0) = x^{0) = ^(0), 
^2(0) = ^4(0) = 0, Pj{0) = 0, 1 < j < 5. We compute the GALIs of an ^5 stable periodic orbit (Figs. [3l^a) 
and (b)) with initial condition x(0) ^ 1.035 for = 5 and an S4U1 unstable periodic orbit (Fig. [3](c)) 
with initial condition x{0) ^ 1.168 for = 7. From Fig. [3] we see again that the behavior of the GALIs 
is well reproduced by Eq. (l22ll for = 5 in the case of the stable orbit, and by Eq. (fMll for ai = 0.088, 
c^i = 0, 2 < z < 5, which are the values obtained by Eq. (fT3ll for the unstable orbit. 
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Fig. 3. The time evolution of (a) GALI2, GALI3, GALI5 and (b) GALIe, GALIg, GALIio for a stable periodic orbit of the 
5dof Hamiltonian system (|26|) . (c) The time evolution of GALI2, GALI3, GALI5 for an S4U1 unstable periodic orbit of the 
same model. Plotted lines correspond to appropriate power laws (|22p in (a) and (b), and exponential laws (|14p in (c). 

4.3. Numerical results - Symplectic maps 

According to the theoretical arguments of Sect. 14. H the GALIs of unstable periodic orbits of maps should 
exhibit the same behavior as in the case of Hamiltonian flows, i.e. they should tend exponentially to zero 
following Eq. (fT4|) . On the other hand, we have argued that the GALIs of stable periodic orbits should 
remain constant, according to Eq. (|23]) . having a different behavior with respect to Hamiltonian systems. 
To verify these predictions, we now proceed to study some periodic orbits in a 2D and a 4D symplectic 
map. 



(27) 



4.3.1. 2D Henon map 

First we consider the 2D Henon map ( Henonl . 

X — xco^{2tiuj) + + x^) sin(27rc<j) 
y = — xsin(27rcL;) + (?/ + x^) cos(27rcL;), 

where is a real positive constant. The phase space of this map for uj — 0.201 is plotted in Fig. [D^a). We 
consider two periodic orbits of period 5 (i.e. after 5 iterations of the map the orbit returns to its initial 
point): an Si stable orbit (blue stars in Fig. ID^a)) with initial condition (x,y) ^ (0.14175, —0.10366), and 
an Ui unstable one (red crosses in Fig. IHa)) with initial condition (x,y) ^ (0.0622148475,0.1477550294). 
Fig. m^b) shows that the GALI2 of the stable periodic orbit oscillates around a constant positive value, 
in accordance to Eq. (|23]) . We have also verified that the GALI2 of the unstable periodic orbit decays 
exponentially to zero following Eq. (IT4]) with ai — 0.0039. 




-0.2 -0.1 



Log(?) 



Fig. 4. (a) The phase space of the 2D Henon map (|27)) for uj = 0.201. The points of two periodic orbits of period 5, a stable 
(blue stars) and an unstable one (red crosses), are also plotted. The time evolution of GALI2 of the stable orbit is plotted in 
(b). Two particular points of the unstable periodic orbit discussed in Sect. [5] are marked by letters A and B. 
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4.3.2. 4D standard map 

Let us now consider the 4D symplectic map Kantz fc Grassberger , 1988l | 



X — X']^ ~\~ X<2 

4 = ^2 + ff sin(2^xi) - |p sin[2^(x3 - xi)] , ^ 
^4 + sin(27rx3) - ^ sin[27r(xi - X3)] 



X3 — 
x^ — 



1), 



(28) 



X4 + ^ sin(27rx3 



|psin[27r(xi 



which consists of two coupled standard maps, with real parameters Ki, K2 and 

In Fig. Owe plot the evolution of GALIs for an 5^2 stable periodic orbit of period 7 with initial condition 
(xi, X2, X3, X4) ^ (0.23666,0,0.23666,0) for Ki ^ K2 ^ 0.9 and /3 = 0.05. Like in the case of the 2D map 
(j2Z|), GALI2, GALI3 and GALI4 remain constant, oscillating around non-zero values, in accordance to 
Eq. (ESj). 




Log(0 



Fig. 5. The time evolution of GALI2 (red curve), GALI3 (green curve) and GALI4 (blue curve) for a stable periodic orbit of 
period 7 of the 4D map (|28)) . 



5. Dynamics in the Neighborhood of Periodic Orbits 

We now turn our attention to the dynamics in the vicinity of periodic orbits, studying initially the neighbor- 
hood of stable periodic orbits in Hamiltonian systems. As a first example we consider the 2dof Henon-Heiles 
system (|24|) . and in particular the stable periodic orbit of period 5 studied in Sect. 14.2.11 In Fig. [U^b) we 
have seen that GALI2 oc GALI3 oc t~^ and GALI4 oc t~^ in accordance to Eq. (f22l) . We expect that 
small perturbations of this trajectory will lead to regular motion on 2-dimensional tori surrounding the 
periodic orbit. For this kind of motion, Eq. ([7]) predicts GALI2 oc constant, GALI3 oc and GALI4 oc t~^. 
Thus, only for GALI2 a different evolution between the periodic orbit and its neighborhood is expected. 
This is actually true, as we see in Fig. El^a) where the time evolution of GALIj^^, A: = 2, 3, 4 is plotted for the 
stable periodic orbit (red curves) and two nearby orbits whose initial conditions result from a Ay = 0.00793 
(green curves) and Ay = 0.02793 (blue curves) perturbation. The GALI2 of neighboring orbits initially 
follows a GALI2 oc t~^ evolution, similar to the periodic orbit, but later on stabilizes to a non-zero value 
as Eq. ([7j) predicts. From Fig. [6]^a) we see that the closer the orbit is to the periodic trajectory the longer 
the initial phase of GALI2 oc t~^ lasts, and the smaller is the final non-zero value to which the index tends. 

Let us now perform a more global study of the dynamics of the Henon-Heiles system. First, we consider 
orbits whose initial conditions he on the py ^ Q line of the PSS of Fig. [I](a). In particular, we use 7000 
equally spaced initial conditions on this line and compute their GALI2 values, using for each of them the 
same set of initial (random and orthonormal) deviation vectors. In Fig. El^b) we plot the GALI2 values 
at t = 10^ as a function of y. The regions where GALI2 has large values (> 10~^) correspond to regular 
motion on 2-dimensional tori. Regions where GALI2 has very small values (< 10~^^) correspond to chaotic 
or unstable periodic orbits, while domains with intermediate values (10~^ ^ GALI2 ^ 10~^^), correspond 
to sticky, chaotic orbits. We also distinguish narrow regions where GALI2 decreases abruptly to values 
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Fig. 6. (a) The time evolution of GALI2, GALI3 and GALI4 for three orbits of the Henon-Heiles system ()24|) : the stable 
periodic orbit of period 5 studied in Sect. 14.2.11 (red curves) and two nearby orbits whose initial conditions result from a 
/\y — 0.00793 (green curves) and /\y — 0.02793 (blue curves) perturbations of the periodic orbit. Note that curves of GALI3 
and GALI4 overlap each other, (b) The GALI2 values at t = 10^ for orbits with initial conditions on the Py — ^ line of the PSS 
of Fig. [TJa), as a function of the y coordinate of the initial condition, (c) Regions of different GALI2 values on the {H2,Py) 
plane of the Henon-Heiles system ([24]). Each point corresponds to an orbit in the neighborhood of a family of periodic orbits 
(white curve) and is colored according to the log(GALl2) value computed at t = 10^. The black filled circle denotes the stable 
periodic orbit of Fig. [TJb). 

10~^ ^ GALI2 ^ 10~^. These correspond to domains of regular motion around the main stable periodic 
orbits of the system, as e.g. in the vicinity of y ^ 0.3 which corresponds to the stable periodic orbit in 
the center of the main island of stability in the PSS of Fig. [l](a). This behavior appears because GALI2 at 
stable periodic orbits decays following a power law and reaches values smaller than the ones obtained 
for the neighboring regular orbits, where GALI2 tends to constant non-zero values, as we have seen in 
Fig. Ha). 

This information can be directly used to identify the location of stable periodic orbits. In Fig. [6](c) we 
show a color plot of the parametric space {H2^Py) of the Henon-Heiles system (l24ll . Each point corresponds 
to an initial condition and is colored according to its log(GALl2) value computed at t = 10"^. Chaotic 
orbits are characterized by very small GALI2 values and are located in the purple colored domains. The 
deep orange colored "strip" corresponds to the vicinity of a family of stable periodic orbits (this family 
is denoted by a white curve) for which GALI2 attains smaller (but not too small) values with respect to 
the surrounding light orange colored region, where regular motion on 2-dimensional tori takes place. We 
note that, as H2 increases, the periodic orbit changes its stability and becomes unstable for H2 ^ 0.146. 
The point H2 — 0.125, Py = 0.14979, denoted by a black filled circle in Fig. [6](c), corresponds to the stable 
periodic orbit of Fig. [l](b). 

The GALIs of chaotic orbits in the vicinity of unstable periodic orbits can exhibit a remarkable oscil- 
latory behavior. Such an example is shown in Fig. [3 for a chaotic orbit of the 2D map (f27l) with initial 
condition {x,y) = (0.06221484498946357,0.14775502681732178) (point denoted by '0' in Fig.|71(b)), which 
is located very close to the unstable periodic orbit of period 5 discussed in Sect. 14.3.11 (point A in Fig. [41(a) 
and Fig. [3(b)). In Fig. [3(a) we see that the GALI2 of this orbit decreases exponentially, reaching very small 
values (GALI2 ^ 10"-^^), since the two initially orthonormal deviation vectors tend to align (Fig. [3(b)) due 
to the chaotic nature of the orbit. 

The evolution of these vectors is strongly influenced by the stable and unstable manifolds of the nearby 
unstable periodic orbit. In particular, as the chaotic orbit moves away from point A along a direction parallel 
to the unstable manifold (green curve in Fig. [3(b)), both deviation vectors are stretched in this direction, 
and shrunk in the direction of the stable manifold (blue curve in Fig. [3(b)). So, after a few hundreds of 
iterations, while the orbit remains in the proximity of point A (note the tiny intervals in both axes of 
Fig. [3(b)), the evolved unit deviation vectors become almost identical, and consequently GALI2 decreases 
significantly. 

Nevertheless, the angle between the two vectors does not vanish, and starts to grow again when the 
orbit approaches point B of Fig. [3(c), which is the next consequent of the unstable periodic orbit (see also 
Fig. [4](a)). The chaotic orbit approaches point B moving parallel to the stable manifold of point B (blue 
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curve in Fig. El^c)). Now the deviation vectors start to shrink along this manifold, while they expand along 
the direction of the unstable manifold of point B (green curve in Fig. El^c)). This leads to a significant 
increase of the angle between the two unit vectors, as we see in Fig. El^c), and consequently to an increase 
of the GALI2 values (Fig. El^a)). 




Fig. 7. (a) The time evolution of GALI2 of a chaotic orbit of the 2D map (|27|) , with initial condition close to the unstable 
periodic orbit discussed in Sect. 14.3.11 The blue curve shows the y coordinate of the orbit in arbitrary units. Consequents 
of this orbit and of two unit deviation vectors from it in the neighborhood of points A and B of the unstable periodic orbit 
of Fig. |4{a), are respectively plotted in (b) and (c). In (b) and (c) the stable and unstable manifolds of points A and B are 
respectively plotted, while the points of the chaotic orbit are labeled according to their iteration number. 



This oscillatory behavior is repeated as the chaotic orbit visits all consequents of the unstable periodic 
orbit, and is clearly seen in Fig. [Tl^a) where the y coordinate of the chaotic orbit is plotted in arbitrary 
units (blue curve) together with the GALI2 values. The horizontal segments of this curve correspond to 
the time intervals that the orbit spends close to the fixed points of the unstable periodic orbit. During the 
first part of these intervals the chaotic orbit approaches a fixed point, the two deviation vectors become 
different and GALI2 increases, while afterwards, the chaotic orbit moves away from the fixed point, whence 
the two deviation vectors tend to align, and GALI2 decreases. GALI2 reaches its lowest values during the 
short transition intervals between the neighborhoods of two successive points of the unstable periodic orbit, 
which correspond to the short connecting segments between the plateaus of the blue curve in Fig. [Tl^a). 
These oscillations of GALI2 can last for quite long time intervals, but eventually the chaotic orbit will 
escape from the strong influence of the homoclinic tangle of the unstable periodic orbit and GALI2 will 
rapidly tend to zero. It is worth mentioning that abrupt cha nges in the valu es of SALI (which practically is 
GALI2) by many orders of magnitude was also reported in Voyatzis . 2008l | for chaotic orbits of planetary 
systems. 

Up to now, we have described in detail these oscillations of GALI2 in the case of the 2D map (|27|) 
because they can be easily explained, while the deviation vectors themselves can be visualized in the 2- 
dimensional phase space of the map. Interestingly, this remarkable behavior occurs in higher dimensional 
systems as well. In Fig. [8] we show two such examples. In particular, we consider a chaotic orbit of the 2dof 
Hamiltonian system (|24|) . whose initial condition is located close to an unstable periodic orbit of period 7 
with initial condition (x, y.Px^Py) - (0, 0.1282112414, 0.4847338571, 0) (Fig. El^a)), and an orbit of the 3dof 
system (|25]) whose initial condition is near the SiUi periodic orbit presented in Sect. 14.2.21 (Fig. [8]^b)). In 
both panels of Fig. [8] we observe an oscillatory behavior of GALI2, similar to the one shown in Fig. El^a). 
We also point out that in both cases all other GALIs show similar oscillatory behaviors. 



6. Connection Between the Dynamics of Flows and Maps 

In Sect. El we discussed the dynamical equivalence between A^dof Hamiltonian systems and 2(A^ — 1)D 
maps, as the latter can be interpreted as appropriate PSSs of the former. We have also seen that GALIs 
behave differently for flows and maps. In particular, as was shown in Sect. (H they remain constant for 
stable periodic orbits of maps (see Eq. (l23]l ) and decrease to zero for flows, according to Eq. (l22]l . 
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Fig. 8. Plots similar to Fig. [7ja) for orbits of (a) the 2dof Henon-Heiles system (|24|) . and (b) the 3dof Hamiltonian system 
(|25|) . Blue curves show in arbitrary units the y coordinate of the studied orbits on (a) the PSS x = 0, ^ of system (|24p , 
and (b) the PSS z = 0, > of system (1251) . 



The fact that maps can be considered as PSS of flows, however, is the key to understanding this 
difference. So, computing the restriction of the GALIs on the PSS of a Hamiltonian system, or more 
generally on spaces perpendicular to the flow, should lead to behaviors of the indices similar to the ones 
obtained for maps. Actually this approach h as already been s uccessfully applied to other chaos indicators 

by only considering the 



related to the evolution of deviation vectors Fouchard et al.. 2002: Barri* 



' applied t 
ril 12005^. 



components of these vectors which are orthogonal to the flow. 

Using deviation vectors orthogonal to the flow, we indeed obtain the same GALI behavior for stable 
periodic orbits of flows and maps. Now, for stable periodic orbits of flows the GALIs of these vectors remain 
constant, as we see from Figs. El^a) and (b) where these GALIs are plotted for the stable periodic orbits of 
Figs. [U^b) and[2]^a), respectively. These behaviors differ, however, from the ones shown in Figs. [Hb) and 
El^a) where the GALIs of the usual deviation vectors were computed. We note that when vectors orthogonal 
to the flow are used, GALl2Ar of an A^dof Hamiltonian system is by definition equal to zero, because the 
2N projected vectors are linearly dependent on a (2A^ — l)-dimensional space. For this reason GALI4 and 
GALIe are not displayed in Figs. [9](a) and (b) respectively. 
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Fig. 9. The time evolution of (a) GALI/., /c = 2, 3 for the stable periodic orbit of the 2dof system ([24]) presented in Fig.[TJb), 
and (b) GALI/^, /c = 2, 3,4, 5 for the stable periodic orbit of the 3dof system ([25)) presented in Fig. [2l^a), when the orthogonal 
to the flow components of the deviation vectors are used. 



7. Summary 

In this paper, we have explored in more detail the properties of the GALI method by using it to study 
the local dynamics of periodic solutions of conservative dynamical systems. To this end, we have: a) theo- 
retically predicted and numerically verified the behavior of the method for periodic orbits, b) summarized 
the expected behaviors of the indices and c) clarified the connection between the behavior of GALIs for 
dynamical systems of continuous (Hamiltonian flows) and discrete (symplectic maps) time. 
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More specifically, we showed that for stable periodic orbits, GALIs tend to zero following particular 
power laws for Hamiltonian flows (Eq. (|22]) ). while they fluctuate around non-zero values for symplectic 
maps (Eq. (123]) ). In addition, the GALIs of unstable periodic orbits tend exponentially to zero, both for 
flows and maps (Eq. (|8j)). 

Finally, we examined the usefulness of the indices in helping us better understand the dynamics in the 
vicinity of periodic solutions of such systems. We explained how, the fact that GALIs attain larger values 
near stable periodic orbits than on the periodic orbits themselves, can be used to identify the location 
of these orbits. We also observed a remarkable oscillatory behavior of the GALIs associated with the 
dynamics close to unstable periodic orbits and explained it in terms of the stable and unstable manifolds 
of the periodic orbit, showing how the influence of these manifolds can lead to large variations of the GALI 
values by many orders of magnitude. 
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